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We study the growth of small-scale inhomogeneities of the density of particles floating 
in weakly nonlinear, small-amplitude, surface waves. Despite the amplitude smallncss, 
the accumulated effect of the long-time evolution may produce strongly inhomogeneous 
distribution of the floaters: density fluctuations grow exponentially with a small but finite 
exponent. We show that the exponent is of sixth or higher order in wave amplitude. 
As a result, the inhomogeneities do not form within typical time-scales of the natural 
environment. We conclude that turbulence of surface waves is weakly compressible and 
alone it cannot be a realistic mechanism of the clustering of matter on liquid surfaces. 



1. Introduction 

Clustering of matter on the surface of lakes and pools and of oil slicks and seaweeds 
on the sea surface is well-known empirically while there is no theory that describes it. 
Since surface flows are compressible even for incompressible fluids, such theory should 
be based on the general description of the development of density inhomogeneities in 
a compressible flow. Important characteristics of the formation of small-scale inhomo- 
geneities is minus the sum of the Lyapunov exponents of the flow, A. It gives the 
asymptotic logarithmic growth rate of the density on fluid particle trajectories at large 
times. The rate A is minus the average value of the velocity divergence seen by a fluid 
particle, and it is always non-negative because contracting regions with negative di 
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vergen c e have more particles and hence larger statistical weight 
( 200lh ; iRuellej (|l996l . Il997i Il999h . We note that A is also the production rate of the 
Gibbs entropy, so that the condition A ^ can be regarded as an analogue of the sec - 



ond law of thermodynamic s for t he dissipative dynamics 
Falkovich fc Fouxo rJ (12003 



see 



Ruellej (Il996l . ll997l[l999h : 



2004). For a generic flow one has A > and the asymp- 



totic de nsity bec o mes a singular measure, the so-called Sinai-Ruelle-Bowen measure, 
see e.g. Dorfmar] (1999). Fo r floaters t his means that they fo r m a multi-fracta l struc - 
ture on the su r face, see Yu et al. ( 1991); iFalkovich et al. d200l|); iBalkovskv et al ( 2001 ); 
Ruellej (fl999l); I Falkovich fc Fouxonl d2003l 120041); iBec et all d2004l); iBalk et all (120041); 
Eckhardt fc Schumacher! (l200lh for theory andlRamshankar et all (ll990l);ISommerer fc Ot 
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1993); So mmererl (11 996): Na meson et all (119961) ;ISchroder et all (|l996f ); ICressman fc Goldbur 
2003); iDenissenko et all (|2006l ); iBandi et all (|2006f ) for experiments. This structure is 
the attractor of the two-dimensional dissipative dynamics obeyed by the part i cles o n 
the s urface , and it evolves constantly for time-dependent flows, see e.g. iDorfmanl ( 19991 ): 
Ottl (|2002l ). The Kaplan- Yorke dimension of the attractor, Dky = 1 + Ai/| A + Ai|, is 
between o ne and two, assuming that the principal Lyapunov exponent Ai is positive, see 
Ottl (l2002t l. 
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Surface flows are generic compressible flows for which the Eulerian compressibility, 
measur ed by the dimensionless ratio C — ((diVi) 2 ) / ((djVi) 2 ) , is of order one, cf. lBoffetta et 
(|2006al ). Here angular brackets stand for the spatial average, v is the floaters velocity 
field, and C changes from zero for incompressible flow to one for potential flow. One ex- 
pects A ~ Ai from C ~ 1, so that the deviation of Dky from the surface dimension 2 is 
also of order one, 2 — Dky ~ 1- The expectation holds true for the flow on the surface of 
three-dimensional turbulence. Pe r forming numerical simulations with the full 3d Navier- 
Stokes equations, Boff etta et al. (12006 q|) found C ~ 5, Dky ~ 1-15 and observe d 
strong clustering on the s urface , see also ICressman et ~al. ( 2004 ); Boffetta et al. ( 200661 ); 
Eckhardt fc Schumacher! ( 2001 ). However, underwater turbulence is relatively rare in the 



natural environment (due to stable stratification), and it is important to consider other 
cases of the surface flows, of which the small-amplitude surface waves are probably the 
most wide-spread. Despite the amplitude smallness, a small-but-finite A produces a large 
effect over time-scales of order 1/A and larger. Thus to estimate the role of surface waves 
in the formation of the floaters inhomogeneities in the natural environment, one needs 
to know how small A is. In this Article we show that A is of sixth or higher order in 
wave amplitude. Let us note that for surface waves the degree of compressibility C is due 
to linear waves, which produce potential flow with C = 1. Thus one could expect that 
the estimates A ~ Ai and 2 — Dky ~ 1 would hold for surface waves like they do for 
the underwater turbulence. We show that in fact, under some natural non-degeneracy 
assumptions described in the Discussion, for surface waves A <C Ai. 

The calculation of A for random waves in different situations was considered in iBalk et al 



( 2004h: IVucelia et all (|2006h . For surfa ce waves, in the li near order in the wave amplitude, 
the particles move periodically, see e.g. lBatchelorl(|2002l ). hence there is no net clustering. 
Thus analysis of clustering of the floaters demands the account of the nonlinear effects. 
Balk et al. ( 20041 ) assumed a linear relation between the velocity field of the floaters and 
the wave amplitudes, and considered a Gaussian ensemble of non-interacting waves. The 
nonlinearity in this case comes from expressing the Lagrangian objects in terms of the 
Eulerian ones. It was shown that A vanishes in the fourth order in the wave amplitude 
for longitudinal waves, whose dispersion relation does not allow the same frequency for 
two different wave- vectors (e.g. sound, gravity, capillary waves). On the other hand, the 
lowest order non-vanishing contribution into Ai was shown to be always of the fourth 
order in the wave amplitude. Under th e same assumption o f the linear relation between 
the velocity and the wave amplitudes, I Vucelja et M (|2006h demonstrated that the ac- 
count of the wave interactions does not change the conclusion of IBalk et all (|2004l ) on 
the vanishing of A in the fourth order in the wave amplitude. 

The results above are inconclusive as far as surface waves are concerned, for which the 
relation between the velocity and the amplitudes is nonlinear due to a small but finite 
curvature of the surface, see IZakharovl (jl966l . 1968 ). For this case a separate calcula- 



tion of A is needed. Here we provide such a calculation. We consider weakly nonlinear 
surface waves and show that neither the wave interactions, nor the nonlinear relation 
between the velocity and the amplitudes, create non-zero sum of the Lyapunov expo- 
nents up to the fourth order in the wave amplitudes. The main tool of our analysis is a 
recently derived Green-Kubo typ e formula for the sum of the Lyapunov exponents, see 
iFalkovich k Fouxonl (|2004i l2003f ). This formula expresses A in terms of the correlations 
of the flow divergence in the particle frame. It describes the interplay between the par- 
ticle motion and the local flow compressions in accumulating density inhomogeneities 
which become pronounced as a result of the long-time evolution. The slowness of the 
particle drift from its initial position allows us to express the correlations in terms of the 
Eulerian correlation functions of the velocity, which we further evaluate by a lengthy, 
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but straightforward calculation, cf. Vucelia et all ( 2006h . The following text is organized 
as follows. In the next Section we introduce the expression for A in terms of the Eule- 
rian correlations of the surface flow velocity, valid up to the fourth order in the wave 
amplitude. To perform the calculation one needs to express the velocity in terms of the 
normal coordinates. This is dealt with in Section [3] The calculation of the various terms 
occurring in A is performed in Section [U Discussion finishes the article. 



2. The sum of the Lyapunov exponents 

Behavior of the density nit, z) in a velocity f ield v(t, x) is governed by the continuity 
equation dtn + V- (to) = 0, see lBatchelorl ( 20021 ). Introducing Lagrangian particle trajec- 
tories by the equation dtX(t, x) = v[t, X(t, x)] with the initial condition X(0, x) — x, we 
may write the solution for the density as n[t, X(t, x)} = n(0, x) exp[— J Q w[t' , X(t' , x)]dt']. 
Here we introduced w(t, x) = V • v(t, x). We characterize the growth of spatial inhomo- 
geneities by the asymptotic logarithmic growth rate A at large times, defined by 



A 



lim — In 

t — >oo t 



i[t,X{t,x)] 
n(0, x) 



lim — 

f^oo t 



w[t',X(t',x)]dt'. 



(2.1) 



The above limit is well-defined. IDorfmanI (| 19991 ). and it gives minus the sum of the 
Lyap unov exponents of the flow v(t,x). It was shown in iFalkovich &: Fouxonl (|2004 , 
20031 ) that if v(t,x) is a random, spatially homogeneous, stationary flow, then one has 



A = 



dt(w(0,x)w[t,X(t,x)]) . 



(2.2) 



We shall apply the above formula to the case where v(t, x) is the two-dimensional velocity 
field governing the motion of the floaters in a (qu asi-)stationary ensemb le of weakly 
nonlinear surface waves sustained by some forcing, see lZakharov et al. ( 1992 ). We first use 
the amplitude smallness to express the Lagrangian correlation function in (12.21) in terms 



of the velocity correlation functions given in the Eulerian frame. We follow [Vuceli a et~al 



(200G) who considered (|2.2p in the case of arbitrary low-amplitude waves. For dispersion 
relation fife, considering packets with both the wavenumber and the width of order k, 
the correlation time of w can be estimated as fi^ 1 and the correlation length as fc -1 . The 
particle deviation from the initial position, X(t, x) — x, during the period, t ~ j is 
e = kv/VLf, <C 1 times smaller than which allows to expand (|2.2p near x. Performing 
the expansion to order e 4 we find 



A k A2 + A3 + A4 , 

/>OC 

A 2 = / dt(w{0)w(t)) 
Jo 

A 4 E 





dt I dti 
i) Jo Jo 



dt 

Jo 



1 /„% . . ( dw(t) dv a (ti) d 2 w(t) , . 

dt, { w(o)v,(t 2 ) f u ;^ + Q^-Mh 



(2.3) 
(2.4) 

(2.5) 



Here and below we suppress the spatial coordi nate x over which th e averaging is per- 
formed. The expansion above was introduced in Vucelia et al. ( 2006f) . We mention that 
all contributions A, are of fourth or higher order in wave amplitude, see below. To use the 
above formula to find A to order e 4 , we need to establish the expression for the surface 
flow v to order e 3 . 
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3. The velocity field of the floaters 

The velocity field that governs the evolution of the floaters coordinates r — (x,y) in 
the horizontal plane, has the following form 

/ .\ (d<j>{r,z,t) d<j)(r,z,t) . \ 
V ^ r ' ' = I dx ^ = ^ ^ ' dy ^ = r| ^ r, " ) ' ' ^ 

where v(r , t) is t he su rface elevation and <fi(r,z,t) is the velocity potential, v = V<^>. 
Zakharov (1966, Il968f) showed that the system of weakly interacting surface waves 



is a Hamiltonian system with canonically conjugate coordinates rj(r,t) and ip(r,t) 
cf>(r, rj(r, t), t). Up to the order e 3 , the surface flow of the floaters is given by 



fciVi - i / T^e^+^lhlk^tm (3.2) 



-\ J^e^+^+^lhl (\h\kS + \h\kg - 2^ + ^)^2773 , (3.3) 

cf. IZakharovl (Jl968) formula (1.8). We introduced shorthand notations rji(t) = n(ki,t), 
4>i(t) = ip(ki,t) and dkiji.., — dkidkjdki.... Note that the velocity field on the surface 
v(r,t) is neither potential, nor solenoidal. We now return to the expression for the sum 
of the Lyapunov exponents (12. 3|) . 



4. Calculation of the sum of the Lyapunov exponents 

In this Section we provide the calculation of A to the fourth order in wave amplitude, 
based on the calculation of the different contributions A;, see (|2.3|) - (12.5I). Some part s 
of this analys i s dea l with the same objects as those considered in iBalk et all (|2004f ); 
Vucelia et al. ( 20061 ). however our analysis is different and it uses specific properties of 
surface waves. Below we prov ide a more detailed c alculation, as compared to the short 
articles IBalk et all (|2004) and I Vucelia et all 1)20061 ) . 



4.1. Contribution of the four-point correlation functions 

We first consider the contribution A4 to A. To the fourth order in wave amplitude we 
may assume Gaussian non-interacting waves and use Wick's theorem to decouple the av- 
erages. Employing identities like (v a (ti)d a dpw(i)) = —({dpv a (ti)) (d a w(t))) 7 that follow 
by integration by parts, one finds 



A 4 = - / dt dh dt 2 
'0 J a Jo 



«<(°)^)M<0»fe)> + (»(0)^ll)(»( t )^H) 



^* a )^W(tiMo)) 



(4.1) 



Here we do not assume isotropy of the waves. Isotropy would make terms like iv a {t\)w{t%)) 
vanish. In the considered order, v = Vip is a potential field and spectral representation 
of the pair-correlation function gives 

f dk f dk 

(V(O)^(t)) = / j^E{k) cos (Q fe t) {v a (0)v fj (t)) = / ——k a kpE{k) cos (fl h t) , 



and similar expressions for other correlation functions in (|4.ip . Note that the potentiality 
of surface waves, holding in the Gaussian approximation, makes the velocity spectrum 
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vanish at k = even if E(k = 0) ^ [cf. iBalk et all (|2004 ): IVucelja et al\ (|2006l )]. We 



find 



dkdqE(k)E(q) 

^O 3 



/ dU 



du 



k q (k ■ q) sin [Slfet] sin [ft g (*i - £2)] 



—/eg (fe • q) cos [fifeii] cos [fiq(t — £2)] — k (k ■ q) cos [fifet] cos [O g (ti — t 2 )] 



k q (k ■ q) sin [Qfe(t - £2)] sin [figti] 



(4.2) 



To calculate the time- integrals we represent the products of trigonometric functions above 
as sums or differences of cosine functions and use 



dt I dt\ I dt2 cos(at + bt\ + ct-i) = — 



ir5(a) ir5(a + b) ir5(a + b + c) 



dt / dt\ I dt2 cos(at — bt\ + H2) : 



b(b + c) be c(b + c) 

n5'(a) w5(a) nS(a — b) 



b 2 



b 2 



b^-c, 



(4.3) 



All terms which are supported only at the zero frequency in the frequency representation 
(5— functions or their derivatives) are also supported only at the zero wavenumber, since 
the dispersion relation of surface waves vanishes at k = only. As a result, due to the 
presence of positive powers of k in (|4.2[) . these terms vanish (similarly to the vanishing of 
the velocity spectrum at k = shown above) . It is then easy to see that the first and the 
fourth terms in A4 (having the same dependence on the wavevectors) cancel each other, 
while the second and the third terms give 



A 4 = 



2^2 



(4.4) 



Since for surface waves the equality of the frequencies of two waves implies the equality of 
their wavelengths, then the above t erms cancel each other, A4 = 0. This reproduces in a 
simple way the result of Balk et al\ ( 2004 ) which says that in the Gaussian approximation 
A vanishes for potential waves having the property that the equality of the frequencies 
implies the equality of the wavelengths. 

4.2. Contribution of the two-point correlation function 

In the next two subsections we calculate those terms in A2 and A3 that can be found 
directly using Wick's theorem. After that we calculate the terms involving the wave 
interactions. We first consider A2 in (|2.3p . Using the explicit form of the surface velocity 
(|3.3[) we get 



2A, 



kfk^dk^dt 



(2^) 4 
dki234dt 



^g^lfclCfcr . k 2 + kl){M0)V2(0)Mt)) 



(2tt) 5 



(^1 (0)772 (0)773 (0)^4 (*)) ^4 1 fci I • k 2 + K + k 2 ■ fc 3 ) 



fell 



2^k\ + kl 



(kx ■ fc 3 + fc 2 • k 3 + k%) 



{2-Kf 



x|/ci||fc 3 |(fci • k 2 + kl){k z ■ k 4 + k\) 



(4.5) 



We note that the third term above, that can be decomposed by Wick's theorem, vanishes 
because it is supported at the zero frequency CO4 imposing k 4 = 0. The last term can also 
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be analyzed using Wick's theorem. Noting that the pair-correlations are supported at 
the zero sum of the involved wavenumbers, we find that the last term equals 



[(M0)Mt))(v2(0)v4(t)) + (M0)v4(t))(m(0)Mt))} 



(27T) 

|fcl||fc 3 |(fcl-fc2 + *a)(fe3-fc4 + *4)- 



(4.6) 



To establish the expressions for the pair-correlation functions we pass to the normal 
coordinates a(k, t): 



n(k,t) 



i(k, t) + a*(-k, t)} , ^(fe, t) = -J ^| [a(k, t) - a*(-k, t)}. (4.7) 



In the Gaussian approximation the pair correlation functions are given by 

(a*(k, t)a(k', 0)) = (2n) 2 5(k - k')n(k) exp [ifl k t] , (a(k, t)a(k' , 0)) = 0, (4.8) 
(i>(k, t)i>{k', 0)} = n k(27T) 2 S(k + k') ^ ^ ( _ , fifet) + n{ _ k) ^ {m _ kt)] ^ 

( v (k, tHk', o)} = fc ( 2 ^) 2 y + fc/ ) [n(fe) exp (_ infet ) + „ ( _ fc ) exp ( in _ fct)] 5 

(^(fe, *)*»(*', 0)) = ( 2 ^) 2(5 ( fc + fc ) [ n(fc) exp f-iQ kt ) _ n (_fe) exp (ifl_ fc *)l . (4.9) 

Using the correlation functions above and noting the vanishing of the terms containing 
8— functions supported at the zero frequency, we can write (|4.6[) as 



dk 12 



2,,, ,. ,'2V2 ( n i k 2 



kf(ki ■ k 2 



[nin_2#(Oi — + n_in2^(fi-i — ^2)] 



--|fex||fc 2 |(fei • fc 2 + fcf)(fei ' fc 2 + fe?) [nin_ 2 £(fil - fi-2) + ri-i^^-i - fi 2 )] 



= 0, 



where we introduced the shorthand notation: n.(±fei, <) = n±i(t) and flifc; = f2±i. We use 
that for surface waves fife is an increasing function of |fe|. We also use that S— functions 
imply Hi = ^2 and k\ — k 2 . We find that A2 can be written as 



A 2 = 



dk\ 2 dt 



k- 2 Jc 2 



(Vh(o)V>2(*)) - 



dk 3 \k\\k\{k\-k 2 + k 2 
(2?r) d 



-(^1(0)772(0)^3^)) 



(4.10) 



The calculation of the above terms demands the account of the interactions. We now 
consider the terms in A3 that can be calculated by the direct use of Wick's theorem. 

4.3. Contribution of the three-point correlation function 
We consider A3 from (12. 3|) . Using the expression (|3.3[) for the velocity, we obtain 



A, 



dk 



12 3;„2 ; „2, 



(2tt) 6 



k(ki(k 2 -k a ) dt dtidn (0)Mt)Mh)) 



dk 



1231 



dt / dt 1 



fci|(fei ■ fe 2 + ki)k£(k 3 ■ k 4 )(MQ)v2(o)Mt)Mti)) 



(2^) 8 Jo 

+k\\k 2 \{k 2 ■ k 3 + k 2 z ){k 2 • fc 4 + fc 3 • k i )(tp 1 {°)M t )v3(t)M t i)) 

+klk 2 2 {k 2 ■ k 4 )\k 3 \(M0)Mt)Mti)V4(h)) • 



(4.11) 
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One can use Wick's theorem for the last three terms. Employing the identity 



dt / dti exp[iujit + iu> 2 t\] 



in[5(uji) - S(uji + ui 2 )] 



W2 



(4.12) 



one finds that part of the obtained terms contain S— functions supported at the zero 
frequency and they vanish because of the vanishing of the integrand there. Let us consider 
the rest of the terms. For the first of the fourth-order terms in (|4.11|) one finds 



dfc 12 f2l7T 



(2?r) 4 4jfci 
|fei|fc 2 (fei • k 2 + k 



(ki-k 2 ){ \k 1 \ 3 (k 1 -k 2 + k 2 2 



n_!n 2 (5(f2_i — f2 2 ) n_ 2 ni5(fl_ 2 — Qi) 



n-in 2 S(ri-i — il 2 ) n^ 2 ni5(fl- 2 — ill) 



= 0. 



where the equality of the frequencies implies the equality of the wavelengths. Analogously, 
for the second Gaussian term one finds 



dk 23 Q 2 tt 
(2tt) 4 4fc7' 



2\2 



\k 2 \ 6 (k 2 -k 3 + ki) 



n_ 3 n 2 5(ri-3 - Q 2 ) n_ 2 n 3 S(£l_ 2 - Sl 3 ) 



k 2 \kl{k 2 ■ fc 3 + k 2 3 ){k 2 ■ k 3 + k\ 



0-3 

n- 3 n 2 S(fl- 3 



n 2 



n 3 

n^ 2 n 3 5(ri- 2 



n 3 



fir 



0. 



Finally, the third Gaussian term contains only 5— functions supported at the zero fre- 
quencies, so it also produces zero. 

4.4. Contribution of the interaction terms 

As a result of the calculation in the previous sections, adding equations (|4. 10|) and (|4.11[) . 
one can write A as a sum over terms which calculation involves the wave interactions. 
These terms have a special property important for the following analysis: all the terms 
contain a field at the zero frequency. Using Fourier representation over the frequency one 
finds 

k\k\dk\ 2 dio 



(2tt) 5 
dki 23 dioi 2 



-(^i(w)V'2(w = 0)) 



(2tt) 8 
dkx^duix 



|fci|fc|(fci • k 2 + k£)(i>i((jj 1 )ri2((jJ2)'ip3(w = 0)> 



23 ,„2 ,2 



(2tt) 9 



k x k 2 (k 2 ■ k 3 ){il)i{uj\)il) 2 (uj 2 )il) 3 (u) 3 )) 



u> 3 



(4.13) 



where we introduced shorthand notation duJiji = du>idu)jdu>i... and used in the last term 
the identity (|4.12[) and the proportionality of the correlation function to 6(lui +uj 2 +u> 3 ). 
We observe that all the terms contain ip(k, ui = 0). We assume that the forcing sustaining 
the stationary wave turbulence vanishes at the zero frequency [note that the first term on 
th e RHS abov e vanis hes i n the Gaussian appro ximation independently of this assumption, 
cf. Balk et al. ( 20041 ) and Vucelja et all (|2006l )] and express ip{k,uj = 0) with the help of 
higher order terms. Substituting the resulting expressions into the correlation functions 
will allow the use of Wick's decomposition. 

4.4.1. The expression for ip(k,uj = 0) 

To calculate A we need to know V'C 8 ? w = 0) to the third or der in wav e ampl it ude. Con- 
sider the dynamical equation on the surface elevation 77, see iBatchelorl (|2002l ) : IZakharov 
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( 1968h . 

drj 



()j . } L z = n]-S/r)V4>[z = r]]. (4.14) 
To order e 2 , the equation reads 

I = tstI* = 0] + = 0] + tst - v? ^ o[z = 0] + 0(7?3) ' (415) 

where </>o and 0i ar e the terms of t he expansion of the potential wi th respect to th e 
surface elevation, see lZakharovl (jl968l ). Using the expressions for 4>i, see lZakharovl (jl968 ). 
performing the Fourier transform over the space and time coordinates, and setting the 
frequency u = 0, we find 

dkidtu 7 2 (/l \ i n f dkiduj , 



0=|A#( fe)CJ = ) + y -^—kt^{ki,u)n(k-ki,-w)-J ■^33-|*||*i|^(fci,w) 
x» ? (fc-fc 1 ,-w) + y ^^fei-Cfe-fei^Cfei^Mfe-fe^-^+O^ 3 ). (4.16) 



It follows that 



^(fe,cj = 0)= / f|fci|-^TT^ ) ^(fci.wMfc-ibi.-wJ + O^ 3 ). (4.17) 



(2tt) 3 

The physical meaning of the above representation is that the zero frequency field arises 
due to the nonlinear interactions only, and in the lowest order it can be represented 
as a result of a single wave scattering. The above formula is sufficient to calculate the 
interaction terms containing the products of three fields. 

4.4.2. Interaction terms containing the products of three fields 

To calculate the interaction terms of the third order we use (I4.17p . Wick's decomposi- 
tion and Fourier transformed version of (|4.9[) . The second term in (|4. 13|) can be written 
with the help of P~T7j) as 

dkx234,du>123 i |,2/, , , i2\{\i. I k 4 ■ k 



(27T) 11 i fel l fe 3( fe i • fe 2 + k 2 ) [\ki\ i^pj (V'(fei,wi)r?(fe2,w 2 )V'(fe4,W3) 

dkuduii , ,,, , 2 /, i2\fn i fci • (fci + fc 2 ) 



77(fc 3 - ^4, -ws)> = - j ^^ IfciKfci + fc 2 ) 2 (fc! • fc 2 + fc 2 2 ) (jfcl 



filial 



\ki + k 2 \ 

[m6(u}i + ili) + ri-\8{u)\ - fl-i)] [n 2 S(—uJi + £1 2 ) + rt-2<5(^i + fi-2)] 



4fi 2 |fci 

1 f dk 12 duj! 2 , 2 \ 1 fc 2 -(fei+fc2 

y (2tt) 7 i fc iK fc i + fc z) ( fc i • fc 2 +^2) - 



|fcl +fc 2 | 

x [ni5(wi + Oi) - n_i5(wi - fi_i)] [n_ 2 £(wi + fi-a) - n 2 <5(wi - fi 2 )] , (4-18) 

where to show that the remaining contraction vanishes, one should use that (ip(k, lu = 0)) 
is representable as an integral over (ip(k, t)rj(k , t)) which in the Gaussian approximation 
vanishes by (|4.9p . Using fi_fc = Qk and noting that the terms supported at fife = 
vanish, one may rewrite the above as 



1 f dkv> , , , ,, , ,0 , /., . fei • (fci + feo) 

7 (2^' l|( 1+ 2) ( 1 2+ 2) (J 1 



4 ./ (2;t)' 11 1 " V |fel+fe: 

1 /" dfe 12 , ,. X2/I. t , l^/i, . 1 fc 2 -(fcl + fc 2 ) 



+n_in 2 <5(ri2 - H-i) 



Various phenomena in fluid mechanics 



[ni»_ 2 a(ni - n_ 2 ) + n_!n 2 5(n 2 - = i y ^^i^Kfex + fc 2 ) 2 (fc 1 • fc 2 + fc 2 2 ) 

(fci - fcg) • (fci + fc 2 ) 
\ki + k 2 \ 



riin_ 2 (S(r2i — 57_ 2 ) + n_in 2 (5(f2 2 — f2_i) 



0. 



(4.19) 



where the last term contains k\ — k\ = 0. Let us now consider the last term in (14. 13)) 
that can be written as 

'' ky ' ''' U23 k\klk z ■ (fc 2 - fei)#^ = 0)4>{k 2 ,uj 2 )4>{k 5 ,u 3 )). (4.20) 



^ 3 (2^) 9 
Substituting (j4~T7j) we find 

■dfci 2 34dw 2 34fc 2 fe|fc 3 -(fc 2 - fci) 



Xlplki, LJi)) = in 



co 3 (2tt)^ 

dk 23 duj 3 



|*4 



fc4 • fci 



uJ 3 {2nf 



(fc 2 + fc 3 )'fc 2 'fc 3 • (2fc 2 + fc 3 ) \k 



(i](ki - fc 4 , -W4)-0(fc 2 , W 2 )^(fc 3 , LU 3 ) 

fc 3 -(fc 2 + fc 3 )\ ft 3 



lfc 2 + fc 3 



4i|fc 3 



[n 2 <5(w 3 - fi 2 ) - n -2$(u 3 + ^-2)] [«3^(^3 + O3) + ™-3<5(w 3 - ^-3)] + i7r 

x(fe+fe3) ^ 3 .( 2fe+ %)( lbl -^a^))(j^) 

x [n 2 <5(w 3 - Sl 2 ) + n_ 2 <5(cj 3 + Q_ 2 )] [n 3 <5(u;3 + ^3) - n_ 3 (5(o; 3 - f2_ 3 )] , 



dk 23 dw 3 
w 3 (2^) 8 



(4.21) 



where the term involving the contraction (^(fci — fc4, — Wijipiki, lo^)) corresponding to 
(tp(k,Lu — 0)}, vanishes as was shown in the analysis of the previous term, where we 
omitted the terms supported at the zero frequency. Taking the integral over u> 3 , omitting 
the terms supported at the zero frequency, we find 



dk 



23 



(27T)* 



(fc 2 + fc 3 )^ 2 fc 3 -(2fc 2 + fc 3 ) \k. 



n 2 n- 3 5(fl 2 - Cl- 3 ) 7l_ 2 Tl3<5(f2_2 - 3 ) 



4i|fc 3 



4i|fc 3 



x(fc 2 +fc 3 r^fc3-(2fc 2 +fc 3 ) ifc 2 |- 



fc 3 ■ (fc 2 + fc 3 
|fc 2 + fc 3 | 

rffc23 

(2tt) 8 



fc 2 • (fc 2 + fc 3 ) 



|fc2 + fc3 

n_ 2 n 3 <5(£l 3 - fi_ 2 ) 7i 2 n_ 3 <S(f2 2 - fl- 3 



M\ki 



4i|fc 2 



= 0, 



(4.22) 



where zero is obtained in the same way as with the previous term. The result of this 
subsection is that A can be written as 



A 



1 f k\k\dk 12 



(2tt)4 



2^ 



(V>(fci,wW(fc 2 ,w = 0)>. 



(4.23) 



4.4.3. Interaction term containing the product of two fields 

To calculate the RHS of (|4.23[) we note that it is sufficient to know i/j(ki,u>) at an 
arbitrarily small but finite u) where the forcing is again negligible and one use the dynamic 
equation (|4.14p without the force. The Fourier transform of (|4.15|) taken now at a small 
but finite frequency gives 



iLurj(k,Lu) = \k\ip(k,u}) + 



dkidtoi 7 2 



(2ir) s 



fc 1 i/'i(wi)?7(fc - fci, w - u>x) 



dk\duj 



-\k\\h\M^) 
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ir 1 {k-k 1 ,u>-u 1 ) + J ^^^-(k-kM^Hk-kt^-^+Oir, 3 ). (4.24) 



which produces 

iujr](k,u)) 



\k\ 



dkiduji 
(2tt) 3 



fell 



fei • k 



tf)(ki,cJi)r](k —ki,uj — uj\) + 0(r] 3 ). 



It follows that the (I4.23|) can be written as 



If k\k\dk\i 



du> I 

(2tt) 4_ J 2^\ 
XT](ki - k 3 , w - wi) + 0(?? 3 ) 
fc 3 • fci 



iur)(k u u) 
ip(k 2 ,u = 0) 



dk 3 dcui 
(2tt) 3 



fca 



fc3 ■ fcl 

l*l| 



1 /" k\k\dk\2 f duj\ 
2 



(2tt) 4 



2tt 



l*i| 



ip(ka, uJ3)v( k i ~ fc 3, wi - u 3 ) + 0(?y ) 



^(fe 2 ,w = 0) 



dkzdtoz 
(2tt) 3 

(4.25) 



where we used that w(^(fei, cj)^/)(fe2, w = 0)) oc oj<5(o;) = 0. Finally, using that ^(fe, = 0) 
is by itself of order e 2 , see (14.17^ . taking the integral over ui 3 and omitting the terms 
supported at the zero frequency, we find A to the order e 4 



(fcl + fe 3 ) 4 |fc 3 | - 



8 7 (2tt) 3 

+n_ a »i*(n 1 - J2_ 3 )] ( |fc 3 |-|fci|- 



|fci +fc 3 

(fc 3 - fc x ) • (fcj + fc 3 ) 
Ifcl +fc 3 | 



[n 3 n_ic)(f2_i — 3 ) 

= 0. (4.26) 



The above can be easily verified by noting that the integrand can be reduced to an 
expression proportional to k\ — k\. Thus we obtain that the sum of Lyapunov exponents 
for weakly interacting surface waves is zero in the fourth order in the wave amplitude. 



5. Discussion 

We have shown that the sum of the Lyapunov exponents for surface wave turbulence 
vanishes in the fourth order in wave amplitude. Using the approximate Gaussianity of 
the statistics, it is easy to see that the leading order term in A is of the sixth order in 
wave- amplitude, or higher. Therefore we have derived that A < fifce 6 . For waves with a 
typical period of the order of seconds and with not too small e = 0.1, we find that the 
time-scale l/fi^e 6 is of the order of weeks. Thus, even if there is no degeneracy in the sixth 
order and A ~ fifce 6 , the formation of the inhomogeneities would occur at the time-scale 
of weeks and larger. It is unlikely that a low-amplitude wave turbulence would persist 
for such time. Thus we expect the turbulence of small-amplitude surface waves to have 
negligible effect on the formation of the floaters inhomogeneities in realistic situations. 
Let us stress that this weak compressibility of surface waves holds in the sense of the 
long-time action of the flow on the particles, while the characteristic value of the ratio 
of the surface flow divergence to the curl is of order one. 

Let us consider some estimates for the Lyapunov exponents of the floaters velocity 
field. For the non-int eracting Gaussian surface waves, Ai is non-zero in the fou rth order in 
wave amplitude, s ee Balk et al. ( 2004 ). while A is non-zero in the sixth order, Balk et all 
( 2004 ); Stepanov ( 2006f h The nonlinear wave interactions and the nonlinearity of the 
velocity-amplitude relation add to Ai additional terms of the fourth order in the wave 
amplitude and higher. It is highly implausible that these terms produce exact cancelation 
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of Ai in the fourth order - such cancelation would depend on the precise form of both the 
interactions and the velocity-amplitude relation. Moreover, positive Lya punov exponent 
andLagrangian chaos hold even for rather simple Eulerian flows, see e.g. lBohr fc Hansen 
(|l996l) . so no degeneracy for Ai is expected that would lead to an exact cancelation in 
the fourth order. Therefore we expect that Ai for surface waves is of the fourth order 
in the wave amplitude. Similarly, we expect no exact cancelation of the non- vanishing 
Gaussian terms for A, in the sixth order in the wave amplitude, see also below. Then we 
have the following order of magnitude estimates: Xi oc f^e 4 and A oc fifce 6 . It follows that 
surface wave turbulence is also weakly compressible in the sense of the ratio A/Ai <C 1, 
which is of the second order in the wave amplitude. The Kaplan- Yorke dimension of the 
particles attractor on the surface is close to the space dimension: Dky ~ 2 — A/Ai i=s 2. 
For dimensionless exponents Ai = Xi/Qk and A = X/Qk we find the order of magnitude 

~ ~ 3 /2 

relation A ~ A-/ holding at small wav e ampli t udes. The above estimates are supported 
by numerical simulations performed bv lUmeki ( 1992[ ) in a similar problem with standing 
surface waves. It was shown there that both A ^ Ai and 2 — Dky -C 1 hold. Moreover, 
the numerical values of the dimensionless exponents A, Ai found there, can be easily 

~ ~ 3 /2 

seen to be in agreement with the relation A ~ A-/ . The box-counting dimension of the 
particles attractor on the surface was found very close to Dky , which we expect to hold 
for wave turbulence as well. The detailed calculations of the exact expressions for Ai and 
A are subjects for future work. 

We believe that our conclusion on the weak clustering in surface waves with a small 
amplitude is an important step in identifying possible reasons for the clusters of the 
floaters observed on liquid surfaces ubiquitously. Our results suggest that other mecha- 
nisms need to be explored, such as the wave breaking and Langmuir circula tion, see e.g. 
Thornel lj20f)5lh the interplay of waves and currents, see lVucelja et all (|2006l ). and others. 



We are indebted to G. Falkovich for constant help and useful discussions. We are 
grateful to M. G. Stepanov for helpful discussions. M. V. thanks S. Nazarenko for a 
useful remark. This work was supported by the Israeli Science Foundation. 
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